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Abstract 



In this note, we establish optimal lower and upper Gaussian bounds 
for the density of the solution to a class of stochastic integral equations 
driven by an additive spatially homogeneous Gaussian random field. The 
proof is based on the techniques of the Malliavin calculus and a density 
formula obtained by Nourdin and Viens in [14] ■ Then, the main result 
is applied to the mild solution of a general class of SPDEs driven by a 
Gaussian noise which is white in time and has a spatially homogeneous 
correlation. In particular, this covers the case of the stochastic heat and 
wave equations in R d with d > 1 and d < 3, respectively. The upper and 
lower Gaussian bounds have the same form and are given in terms of the 
variance of the stochastic integral term in the mild form of the equation. 
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1 Introduction and main result 



Let X = {X(t,x), (t,x) G [0, T] x R d } be a zero mean Gaussian process with 
an homogeneous covariance function of the form 

E(X(t,x)X(s,y)) = *(t,s,x-y), 

where we assume that £ i— » $(£, £,0) is continuous, £, 0) > for any £ > 0, 
and $(t,0,x) = for all (t,x) G [0,T] x R d . 

The purpose of this note is to establish optimal Gaussian lower and up- 
per bounds for the probability density of the solution u(t, x) to the following 
stochastic integral equation: 

u(t,x) = X(t,x) + [ [ b(u(s,x-y))T(t- s,dy)ds, (1) 
Jo Js. d 

where (£, x) £ [0, T] x R d and T(£, dy) is a nonnegative measure such that 

sup r(£,M d ) < +oo. (2) 

0<t<T 

Note that if b : R — > R is Lipschitz continuous, there exists a unique solution of 
([1} and the process it is also homogeneous is the space variable. 

In particular, these equations include mild solutions to a large class of SPDEs 
with additive noise of the form 

Lu(t,x) = b(u(t,x))+aW(t,x), (£, x) G [0, T] x R d , (3) 

where L denotes a second order differential operator with constant coefficients, 
the random perturbation W(t, x) stands for a Gaussian noise which is white in 
time and has a spatially homogeneous correlation, and a is constant. We as- 
sume here vanishing initial conditions. Let us denote by \i the spectral measure 
associated to the noise W (for the precise description of the noise, see Section 
H]). The process X(t,x) is in this case 

X(t,x) = a [ [ r(t-8,x-y)W(ds,dy), 

JO JR d 

where L denotes the fundamental solution associated to L and, if we simply 
denote $(£) := $(£,£,0), then 

$(£) = / / \JT(s)(0\ 2 »(d£)ds, 

JO JR d 

This has to be a finite quantity for X(t,x) to be well-defined. Let us point 
out that the main examples of SPDEs of the form © to which our results will 
apply are the stochastic heat equation in any space dimension and the stochastic 
wave equation in R d with d < 3. Indeed, these types of SPDEs have been widely 
studied during the last two decades, see e.g. [21 1 l2l [T2l [51 [20 ) IT9] . 
The main result of the paper is the following: 
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Theorem 1 Assume that T satisfies 0) and b is of class C 1 with a bounded 
derivative. Then, there exists Tq > such that for all (t,x) € (0,Tq) x M. d , 
the random variable u(t,x) solution to fl]) has a density p satisfying, for almost 
every z£l: 

E \u(t,x) — ml f (z - m) 2 1 , , £7|«(t, a;)— ml { (z 
■exp<^- v / }<p(z)< I,'' -exp^ 



C 2 $(*) I Ci$(t) J-*^'- Ci*(t) I C 2 $(t) 

where m = E(u(t, x)) and the constants C\, are positive and only depend on 
b and Tq . 

In order to prove Theorem [TJ we will apply the techniques of the Malliavin 
calculus. Indeed, as it has been also done in [18], we will make use of the 
recent results obtained by Nourdin and Vicns in |14j . In this paper, the authors 
provide sufficient conditions on a one-dimensional Wiener functional (that is, 
on a random variable on the Wiener space) so that its law has a density and it 
is given by an explicit formula (see [21 Theorem 3.1 and Corollary 3.3]). As it 
will be made precise in the proof of Theorem[T] (see Section[3|), the application of 
the above-mentioned result of Nourdin and Viens will require a careful analysis 
of the Malliavin derivative of the solution of equation |T]). 

We should also remark that recently, there have been an increasing inter- 
est in applying the techniques of the Malliavin calculus in order to establish 
Gaussian lower bounds for the probability density of a general class of Wiener 
functionals. In particular, this has been applied to diffusion processes and solu- 
tions to SPDEs, and the article by Kusuoka and Stroock [9] can be considered as 
the starting point of this methodology. Therein, a Gaussian type lower bound 
for the density of a uniformly hypoelliptic diffusion whose drift is a smooth com- 
bination of its diffusion coefficient was established. Later on, Kohatsu-Higa |8J 
got rid of that dependence constraint on the drift and, moreover, took the ideas 
of Kusuoka and Stroock in order to construct a general method to prove that 
the density of a multidimensional functional of the Wiener sheet in [0, T] x M. d 
admits Gaussian lower bounds (see [7]). Then, in the latter paper, Kohatsu- 
Higa dealt with a one-dimensional stochastic heat equation in [0, 1] (for a related 
result, see also [T5] Theorem 3.1]) and his result was also applied by Dalang and 
Nualart [4] in the case of a one-dimensional reduced wave equation. The ideas 
of [7] have also been further developed by Bally [1] to obtain Gaussian lower 
bounds for locally elliptic Ito processes. Eventually, another recent method for 
deriving Gaussian lower bounds for multidimensional Wiener functionals has 
been obtained by Malliavin and Nualart [TU] (see [T5] for a one-dimensional ver- 
sion of this result). This method is based on an exponential moment condition 
on the divergence of a covering vector field associated to the Wiener functional. 

The paper is organised as follows. In Section [2j we will briefly recall the 
main tools of the Malliavin calculus needed in the proof of the main result. In 
particular, we will recall the main points of the method by Nourdin and Viens 
[14] . Section [3] is devoted to prove Theorem [T] Finally, in Section [4] we apply 
Theorem [1] to the solution of a general spatially homogeneous SPDE. 
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2 Preliminaries 



In this section, we will briefly describe the Gaussian setting in which we will 
apply the techniques of the Malliavin calculus as well as the method of Nourdin 
and Viens [13]. For a more complete account on these methodologies, we refer 
the reader to [TB] and [HI Sections 2 and 3], respectively. 

We denote by H the Gaussian space generated by the Gaussian process 
{X(t,x), (t,x) <E [0, T] x M d }; for a complete account on Gaussian Hilbert 
spaces, we refer the reader to [6l Chapter 1]. 

Then, as usual, we denote by D the Malliavin derivative, defined as a closed 
and unbounded operator from L 2 (Q) into L 2 (Tl;7i), whose domain is denoted 
by B 1 ' 2 . Thus, for any random variable F belonging to D 1,2 , its Malliavin 
derivative DF defines an element in L 2 (f2;7i). 

Another important operator in the Malliavin calculus theory that plays an 
important role in [Hj is the generator of the Orstein-Uhlenbeck semigroup (see 
[TrU Section 1.4]). It is usually denoted by L and can be defined by means of its 
Wiener chaos expansion: 

oo 
n=0 

where J„ denotes the projection onto the nth Wiener chaos. 

In P3], the authors consider a random variable F £ D 1,2 with mean zero 
and define the following function in R: 

g F {z) := E[{DF, -DL^F^F = z], 

where L _1 denotes the pseudo-inverse of the generator of the Orstein-Uhlenbeck 
semigroup L. We observe that, by [131 Proposition 3.9], it holds g(z) > on 
the support of F. Then, Nourdin and Viens prove the following result (see [T4l 
Theorem 3.1 and Corollary 3.3]): 

Theorem 2 Assume that there exists a positive constant ci such that 

gF{F) > ci, a.s. 

Then, the law of F has a density p whose support is K and satisfies, almost 
everywhere in M; 

p( z ) = ^— tt ex P ( - / —^dy ) ■ 

2gF{z) V Jo 9F{y) J 

As stated in [141 Corollary 3.5], an immediate consequence of the above theorem 
is that, if one also has that gF(F) < C2, a.s., then the density p satisfies, for 
almost all z£l: 



In order to deal with particular applications of this method, [14, Proposition 3.7] 
establishes an alternative formula for gp(F). Indeed, given a random variable 
F G D 1,2 , one can write DF = $>f(W), where $f is a measurable mapping 
from M. n to Tt, determined (P o ^-almost surely (see [TB], p. 54-55). Then, 
it holds that 



9f(F) 



'E 



(G> F (W),$ F (e~ e W+ Vl - e- 2e W')) u \F dS, (4) 



where W stands for an independent copy of W such that W and W 1 arc defined 
on the product probability space {Vt x O' , T <S> T' , P x P 1 ) . Eventually, E denotes 
the mathematical expectation with respect to F x P'. 

Let us observe that formula ((¥]) can be still rewritten in the following form: 



9f(F) 



>E 



E'((DF,DF) n )\F 



66. 



where, for any random variable X defined in (SI, J 7 , P), X denotes the shifted 
random variable in SI x SY 



X(u,oj') = X(e- e uj + Vl-e- 2e uj'), u € fi, J e SI. 



Notice that, indeed, X depends on the parameter 8, but we have decided to 
drop its explicit dependence for the sake of simplicity. 



3 Proof of the main result 

This section is devoted to prove Theorem[T] To begin with, we have the following 
result, whose proof is straightforward and omitted. 

Proposition 3 Suppose that T satisfies condition JI|) and b : R —* R is of class 
C 1 with a bounded derivative. Let u(t, x) be the solution to Equation (Qp. Then 
u(t,x) belongs to the space D 1 ' 2 and 

Du(t,x) = X(t,x) + / / b' (u(s, x — y))Du(s, x — y)T(t — s, dy)ds, (5) 

a.s. for all (t,x) G (0, T] x R d . 

Let us remark that the pathwise integral in |5j) takes values in the Hilbert space 
7i and can be defined using a standard procedure (see e.g. [El p. 292]). 

Owing to Theorem [5J Theorem Q] will be a consequence of the following 
proposition. For (t,x) G (0, T] x R d we set F := u(t,x) — E(u(t,x)), so that 
we need to find almost sure upper and lower bounds for the random variable 
9f(F): 

g F {F) = J e- e E E' {{Du{t,x), Du(t,x)) n ) \F 60. (6) 
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Proposition 4 Under the same hypothesis as in Theorem^ there exists Tq > 
such that, for all t £ [0,Tq): 

Ci$(t) < g(F) < C 2 $(t), a.s., 

where C\, Ci are positive constants depending on b and Tq. 

For the proof of this proposition, we need the following lemma. Recall that 
:= $(t,i,0) = E(X(t,x) 2 ). 

Lemma 5 There is a positive constant C such that, for all t > 0: 



sup E[\\Du(r,y)\\ H \F] < <V*(*), (7) 

0<r<t 



y& 

and 

sup sup E E 

8>1 o<T<t 



' (||2M^i)|| W ) \f] < Cy/Wfy, a.s. (8) 



Proof. From §E§ and applying Minkowski inequality, we get: 
\\Du(t,y)\\ H <y/m + W\L [ I \\Du(s,y-z)\\ H T(t-s,dz)ds. 

JO JR d 

As a consequence, we have the following estimate: 



E[\\Du(t,y)\\ H \F] < Vm + Wl 

E [\\Du(s, y - z)\\ H \F] T(t - s, dz)ds. 



oo 

t 



^R d 

Set 

Y t := sup E[\\Du{r,y)\\ H \F]. 

0<r<t 
y&R d 

Then, we have proved that 

Y t < VWj+Wb'W^ [ Y s T(t-s,R d )ds, 
Jo 

and a suitable generalisation of Gronwall lemma (e.g. [3l Lemma 15]) allows us 
to conclude the proof. ■ 

Proof of Proposition [4j It follows the same lines as in the proofs of 
Propositions 4.5 and 5.5 in [TS]. Indeed, one just needs to be careful and try to 
keep the sharpest bounds which appear throughout the proof. Let us be more 
precise about this method, namely, by ([5]) and ©, we have that 

3 

g F {F)=$(t)+^A i (t,x), 

i=l 



G 



where 
Ai(t,x) 

A 2 {t,x) 
A 3 (t,x) 



= E 



b'(u{s,x- y))(X(t,x),Du(s,x - y))nT(t - s,dy)di 



e-»E 



E' 



JR d 



b'(u(s,x- y)) 



x (X(t, x),Du(s, x - y)) n T(t — s, dy)dsj 



F 



e-"E 



E' 



b'(u(s, x — y))b'(u(r, x — z)) 



[0,t] 2 JR 2d 

x {Du(s, x — y), Du(r, x — z)}nT(t — s, dy)T(t — r, dz)drds 

Let us first prove the lower bound. For this, observe that we have 

9f(F) > $(t) - \A x {t,x) + A 2 (t,x) + A 3 (t,x)\. 
Let us bound |j4i(i,a:)|. We will use the following notation: 

*(t) := / T( Sl R d )ds. 



dO. 



(9) 



Notice that, by hypothesis, sup tg [ r i $f(t) < +00 and in fact ^(i) converges to 
as t tends to 0. We can write, applying (JT]) in Lemma[5]and using the notation 



\A x {t,x)\ < ^/W) ll&'ILE 



< VW) W\ 



\\Du(s,x- y)\\ n T(t - s,dy)ds 
'0 JR d 

T(s,R d )ds) sup E[\\Du(s,z)\\ n \F] 



F 



0<s<t 

z&i 



For the analysis of the term ^(^x)!, one can proceed using exactly the same 
arguments (the expectation E 1 and the integral with respect to 9 do not affect 
the final result), but here we will need to apply ([5]) in Lemma [5] Indeed, one 
also obtains that 

\A 2 (t,x)\ KCpW^Mt). 
Eventually, in order to bound (^(i, x)\, we observe that 



\A 3 (t,x)\ < 



/||2 

oc 



xE 



[0,t] 2 JR 2d 



\\Du(s,x-y)\\ n E' (\\Du(r, 
xT(t - s, dy)T(t - r, dz)drdsd6. 



x - z)\\ H 



F 
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for which we apply Cauchy-Schwartz inequality with respect to the conditional 
expectation with respect to F and we use the estimate 



sup sup E 

e>i o<T<t 



E 



' (\\Du(r,' t 



\ndr 



F 



< C*(<), 



a.s., 



whose proof is similar to that of ([5]). In this way, we obtain that 

\A 3 (t, X )\<c\\b>\\i<t>m(t) 2 - 

Putting together the bounds for \Ai(t, x)\, i = 1,2,3, we have that there are 
positive constants k\ and k 2 (only depending on b) such that, for all t S (0, T] 
and x e R d (recall that F = u(t,x) - E(u(t,x))), 

g F {F) > - k 2 ($(t)¥(t) + <t>(t)V{t) 2 ) 

= $(i) (h-k 2 (*(t) + *(t) 2 )) , 

where we remind that #(t) — f Q T(s, R d )ds. Hence, if t is sufficiently small, say 
t < T with T > satisfying fc x - fc 2 (*(T ) + *(T ) 2 ) > 0, then 

> Ci*(*), 

where Ci is a positive constant depending on 6 and To. This proves the lower 
bound in the statement. 

Concerning the upper bound, observe that we have, for all t £ (0, To): 



g P (F)<9(t)+^2\Ai(t^)\ 



i=l 

< C$(t) (1 + #(T ) + *(T ) 2 ) 

< <7 a *(t), 

where C2 is a positive constant depending on 6 and T). This concludes the 
proof. ■ 



4 Application to spatially homogeneous SPDEs 

As we commented in the Introduction, the main examples of equations of the 
form (JTJ) which Theorem [T] can be applied to correspond to the mild formulation 
of the following SPDEs: 

Lu(t,x)=b(u(t,x))+aW(t,x), (t, x) £ [0, T] x R d , (10) 

where L is a second order differential operator with constant coefficients and W 
is a Gaussian noise which is white in time and has some spatially homogeneous 
correlation. More precisely, on some probability space (O,^ 7 , P), we consider 
a family of centered Gaussian random variables W = {W(<p), <p G C Q yo (M. d+1 )}, 
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where Cg°(R d+1 ) denotes the space of infinitely differentiable functions with 
compact support, with the following covariance functional: 

E(W{<p)W(il>))= f f {<p(t)*^ a) (t))(x)A(dx)dt, ^VeC °°(M d+1 ), 

JO JR d 

(11) 

where ip^(t,x) := ip(t, —x) and A is a non-negative and non-negative definite 
tempered measure. For the right-hand side of (fTTj) to define a covariance func- 
tional, it turns out that A has to be symmetric and the Fourier transform of 
a non-negative tempered measure fi (see [331 Chapter VII, Theoreme XVII]). 
Usually, \i is called the spectral measure of the noise W. Let us denote by (Tt)t 
the filtration generated by W, conveniently completed. 

Usual examples of spatial correlations are given by measures of the form 
A(dx) = f(x)dx, where / is a non-negative and continuous function on R d \{0}, 
which is integrable in a neighbourhood of 0. For instance, one can consider 
a Riesz kernel f(x) = \x\~ e , for < e < d, while the space-time white noise 
corresponds to consider / = So; in this latter case, the spectral measure is the 
Lebesgue measure on Mr. 

By definition, a mild solution of (fTTTj) is a .Ft-adapted stochastic process 
{u(t,x), (t,x) 6 [0,T] x R d } satisfying 

u(t,x) = <t I I r(t — s,x — y)W(ds,dy) 
Jo Js. d 

+ [ f b(u(s,x-y))T(t-s,dy)ds, (12) 

JO JR d 

a.s. for all (t, x) € [0, T] x R d , where T denotes the fundamental solution 
associated to L. Without any loose of generality, we can assume that <j = 1. Let 
us point out that the stochastic integral in the right-hand side of equation (|12[) 
takes values in R and is considered as an integral with respect to the cylindrical 
Wiener process associated to W, as described in [T71 Section 3]. Note, however, 
that this stochastic integral can be also defined in the sense of Dalang [3] (see 
also [23]) or even in the more abstract framework of Da Prato and Zabczyk [5]. 
In this setting, the process X is given by 

X(t,x)= [ f r(t-s,x-y)W(ds,dy), 

JO JR d 

for which one verifies that 

E(X(t,x)X(s,y))= r [ e 2 ^ x -y^T(t - r)(0W(s - r)(£) »(d£)dr. 
Jo J& d 

Hence, in particular, 

*(*) =E(X(t,x) 2 ) = [ f \TT(r)(m 2 KdOdr, t G [0,T]. 

JO JR d 

At this point, we observe that X(t, x) is a well-defined Gaussian random variable 
whenever the following condition is satisfied (see e.g. [T7J Lemma 3.2]): 



9 



Hypothesis 6 For all t > 0, T(t) defines a non-negative distribution with rapid 
decrease such that 



( I \FT(t)(0\ 2 ti(dt)dt <oo. 

JO JR d 



Moreover, T is a non-negative measure of the form T(t, dx)dt such that, for all 
T>0, 

sup T(t,R d ) <C T <oo. 

0<t<T 

Thus, Theorem [1] applied to equation (fT2")) reads: 

Theorem 7 Assume that Hypothesis [3] is satisfied and b is of class C 1 with a 
bounded derivative. Then, there exists T$ > such that for all (t, x) £ (0, To) x 
W 1 , the random variable u(t, x) solution to ilty) has a density p satisfying, for 
almost every z £ R: 

E\u(t,x)-m\ { (z-m) 2 ) , . E \u(t, x) - ml f (z 



C 2 $(t) ^1 <7i*(t) J-^"'- Ci$(t) C 2 $(t) 

where m = E(u(t, x)) and the constants Ci, C 2 are positive and only depend on 
b and Tq . 

Example 8 Theorem [7] applies to the case of the stochastic heat equation in 
any space dimension and the stochastic wave equation in M. d with d < 3. Indeed, 
in both examples, Hypothesis is fulfilled if and only if 

/Lt(d£) < +oo 



i + iei 

(see e.g. Section 3]). As a consequence, in the case of the stochastic wave 
equation, Theorem^ exhibits an improvement of J18[ Theorem 5.3], where the 
lower and upper Gaussian bounds were not optimal and a slightly stronger con- 
dition on pL was assumed. 

Remark 9 Independently of the results in the study of the existence and 
smoothness of the density for SPDEs of the form ilO\) ( in particular for stochas- 
tic heat and wave equations) , has already been tackled by several authors. Namely, 
let us mention the works JMUMHIWMlMWj, 
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